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ABSTRACT 

This  paper  presents  the  basic  considerations  behind  sight-time 
stabilization.  The  equations  of  motion  that  relate  the  general  pitch,  roll 
and  yaw  motions  of  a  vehicle  hull  to  the  traverse  and  elevation  angles  of 
a  sight-line  fixed  to  the  hull  are  developed.  A  detailed  mathematical  model 
of  a  feedback  control  system  that  will  accomplish  sight-line  stabilization 
is  developed,  and  a  numerical  optimization  is  performed  on  a  simplified 
version  of  this  detailed  model.  The  optimization  routine  used  is  described 
in  reference  [6]. 
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CHAPTER  1 
INTRO  AUCTION 


The  optimization  of  feedback  control  systems  is  a  topic  that 
has  received  much  attention  during  recent  years.  Analytical  techni¬ 
ques  have  been  developed  for  a  great  variety  of  problems,  but  the 
application  of  these  techniques  to  complex  problems  can  become 
extremely  difficult,  if  not  impossible.  This  difficulty  in  apply¬ 
ing  analytical  techniques  has  been  one  of  the  influential  factors 
on  the  rise  of  numerical  techniques  for  optimization. 

Numerical  techniques  can  be  divided  into  «.wo  categories:  direct 
search  and  indirect  search.  All  indirect  methods  employ  the  use  of 
the  gradient  of  the  performance  surface  in  one  manner  or  another. 

For  systems  with  highly  oscillatory  or  even  discontinous  performance 
surfaces,  which  are  frequently  encountered  in  the  feedback  control 
systems  area,  it  is  very  easy  to  calculate  an  incorrect  gradient, 
thus  affecting  the  convergence  of  the  optimization  technique.  An 
incorrect  gradient  is  shi  in  fig.  1. 

Because  of  the  difficulty  involved  in  calculating  the  gradient 
of  the  performance  surface,  direct  methods  are  more  suited  to  control 
optimization  problems.  There  are  many  fine  direct  methods,  among 
them  are  the  method  of  Hook  and  Jeeves  {4}  and  the  method  of 
Rosenbrock  {9}.  T.  Lange- Nielsen  {6}  has  modified  Rosenb rock's 
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method  and  used  it  to  develop  an  optimization  algorithim  using  CSMP  . 

Details  of  Rosenrrock 1 s  method  and  Lange* Nielsen's  modifications 
can  be  found  in  reference  [6], 


When  using  a  numerical  optimization  routine  one  must  be  careful 
to  Insure  that  the  results  are  correct.  In  optimization  the  objective 
is  to  minimize  a  performance  index  while  a  system  goes  from  one  state 
to  another.  For  a  performance  index  such  as  the  integral  squared 
error  (ISE)  an  analytical  check  can  be  made  of  the  results  by  using 
tables  {7}  that  express  ISE  in  terms  of  known  system  quantities. 


Continuous  Systems  Modeling  Program.  An  IBM  analog  simulation 
that  uses  digital  methods,  (see  reference  {5}) 


This  is  not  to  say  that  the  system  can  be  easily  optimized  by  analy¬ 
tical  methods  if  ISE  is  used.  Another  wav  to  help  guarantee  that 
results  are  correct  is  to  be  familiar  enough  with  the  problem  that 
is  being  optimized  that  you  have  some  idea  of  what  your  results  should 
be.  I'c  is  for  this  reason  that  the  material  in  chapters  two,  three 
and  four  is  presented. 

Care  must  be  taken  when  synthesizing  a  mathematical  model.  An 
inaccurate  model  can  cause  all  efforts  to  have  been  in  vain,  and  the 
results  trill  be  useless.  Automatic  control  problems  include  the 
dynamics  of  mechanical  components  as  well  as  the  characteristics 
of  electrical  and  hydraulic  components.  Sight-line  stabilization  is 
a  problem  that,  requires  some  modeling  in  all  three  of  these  areas  as 
well  as  a  thorough  understanding  of  kinematics,  and  of  course  control 
theory. 

The  problem  of  sight-line  stabilization  will  be  examined  in 
chapters  two.  three  and  four;  and  .  detailed  mathematical  model  will 
be  developed  in  chapter  four.  Sometimes  concessions  have  to  be  made 
to  factors  over  which  there  is  no  control.  Computation  time  on  a 
computer  is  one  such  factor.  It  is  for  this  reason  that  a  simpli¬ 
fied  version  of  the  detailed  model  is  used  in  the  actual  numerical 
optimization.  This  simplified  model  is  developed  in  chapter  five. 


Chapter  2 


SIGHT-LINE  STABILIZATION 

The  concept  behind  sight-line  stabilization  is  a  simple  one. 

Consider  a  sight-line  unit  vector,  es>  (see  fig.  2)  directed  from 

a  point,  0,  to  a  point,  P,  with  its  base  fixed  at  0.  Point  P  is 

fixed  in  inertial  space,  and  0  is  fixed  to  a  base,  B,  as  shewn  in 

fig.  2.  The  control  variables  are  the  traverse  angle,  T,  and  the 

elevation  angle,  E,  defined  as  in  fig.  2.  The  problem  is  to  keep 

e  oriented  from  0  to  P  while  the  hull  undergoes  general  pitch, 
s 

roll,  and  yaw  motion  about  the  z,  x,  and  y  axes  respectively. 

The  notation  will  be  to  use  p  to  denote  pitch  angle,  r  to  denote 
roll  angle,  and  q  to  denote  yaw  angle;  with  superscript  dot  denoting 
the  familiar  derivative  with  respect  to  time  (p,  r,  and  c  are  shown 
in  fig.  2). 

It  is  now  desireable  to  derive  an  expression  for  T  and  E  in 
trrms  of  the  known  quantities.  In  the  future  it  will  be  found  that 
it  ir  useful  to  do  this  in  terms  of  the  x-y-z  reference  frame  pre¬ 
viously  defined.  The  sight-line  unit  veetc-r,  eg,  car.  be  vritten 
in  terms  of  the  x,  y  and  z  axis  unit  vectors;  I,  j,  k  respectively. 

e  *  cos(E) cos{T)» i  +  sin(E)-j  +  cos(F,)sin(r)*k  (1) 


(2) 


B 


k  (eV 


+  w  x  e 


s 


t 


The  superscript:  I  denotes  the  inertial  reference  frame,  the  super¬ 
script  B  denotes  the  base  reference  frame,  and  a  is  the  general 
rotation  vector  for  the  base  given  by  equation  (3) . 

(o  =  ri  +  qJ  +  pk  (3) 

By  expanding  the  expressions  on  the  right  of  equation  (2)  equations 
(4)  and  (5)  are  obtained. 

u  x  e  “  {  c  •cosv'E)sin(T)  -  £*sin(E)  }*i  + 
s 

{  p-cos 'E)cos(T)  -  t»cos(E)8in(T)  }*j  + 

{  r-sin(E)  -  q-cos(E)cos(T)  }*k  (4) 

B 

—  (e  )  »•  {  -  T-cos(E)sin('f)  -  E*sin(E)cos(T)  }*i  + 
ut  S 

{  E-cos(E)  }-j  + 

{  T*cos(E)cos(T)  -  E*sin(E)sin(T)  }*k  (5) 


The  stabilization  requirement  is  given  by  equation  (6): 


<e  J  -  0 


(6) 
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Equation  (2)  therefore  reduces  to  equation  (7): 

B,  _  . 

0  =  —  (e  )  +  <d  x  e  (7) 

dc  s  s 

By  substituting  equations  (4)  and  (5)  into  equation  (7)  and  noting 
that  for  a  vector  quantity  to  be  zero  each  component  must  be  zero, 
the  following  set  of  three  scalar  equations  from  the  i,  j  and  k 
components  respectively  are  obtained. 

0  =  -E*cos (T)sin(E)  -  T* sin (T) cos (E)  - 

p-sin(E)  +  q*sin(T)cos(E)  (8) 

0  «*  E*cos(E)  +  p-cos(T)cos(E)  - 

r*sin(T) nos(E)  (9) 

0  =  -E*sin (T)sin(E)  +  T-cos (T)cos (E)  + 

r*sin(E)  -  q»cos (T)cos(E)  (10) 

By  solving  equation  (9)  the  expression  for  E  is  obtained. 

E  »  r*sin(T)  -  p-ccs(T)  (11) 

Substituting  equation  (11)  into  either  equation  (8)  or  equation 
(10)  gives  the  expression  for  T,  equation  (12): 


T  =  q  -  tan(E){r-cos(T)  +  p*sin(T)} 


(12) 


Equations  (11)  and  (12)  are  the  governing  equations  of  motion. 

Note  that  T  and  E  are  measured  with  respect  to  the  base. 

Now  that  the  equations  of  motion  for  T  and  E  are  known  it  must 

be  determined  how  T  and  E  are  to  he  controlled.  Associated  with 

the  traverse  axis  there  is  an  inertia,  J^,  and  with  the  elevation 

axis  an  inertia,  Jg£.  These  are  represented  in  fig.  3.  rotates 

with  respect  to  the  base  around  the  previously  defined  y  axis,  and 

is  thus  measured  with  respect  to  the  y  axis.  Define  a  new  set 

of  coordinates  x^,,  y T,  z^,  that  are  fixed  to  J.^  such  that  y^.  and 

z_  rotate  with  J_.  J__  is  confined  to  rotate  about  the  z_  axis, 

I  T  ££  T 

and  is  thus  measured  vith  respect  to  the  axis.  It  is  easy  to 
see  that  part  of  J£E  will  be  included  in  J^.  This  component  is 
denoted  by  Jg^.  is  not  fixed,  but  is  a  function  of  E.  The 

movement  of  through  an  angle  T  is  accomplished  by  a  motor,  M^., 
mounted  on  parallel  to  the  y  axis  as  shown  in  fig.  3b.  The  torque 
is  transferred  through  a  gear  ratio  of  NT>  Similarly  J£  Is  driven 
by  a  motor,  Mg,  mounted  on  JT  parallel  to  the  z T  axis  with  an 
associated  gear  ratio  Ng. 

It  is  desireable  to  return  to  tht  governing  equations  of  motion 
and  examine  the  significance  of  their  terms  more  carefully.  They 
are  repeated  below  for  convenience. 


E  *  r*sin(T)  -  p-cos(T) 


(ID 


T  ■  q  -  tan(E){r *cos(T)  +  p*sia(T)} 


(12) 


The  pitch  and  roll  velocity  vectors  can  be  resolved  through  the 

traverse  angle  T  to  fora  two  new  vectors,  r*cos(T)  +  p‘Sin(T)  and 

p»cos(T)  -  r*sin(T),  which  together  with  the  original  yaw  vector, 

q,  fora  equivalent  components  of  hull  motion  (see  fig.  4).  It  is 

easy  to  see  that  r-cos(T)  -!  p*sin(T)  is  in  the  previously  defined 

x^,  direction,  and  p-cos(T)  -  r*sin(T)  is  in  the  previously  defined 

z„  direction.  One  of  the  stabilization  requirer*:  its  is  to  drive 

E  about  the  z^  axis  such  that  eg  remains  oriented  in  inertial  space. 

From  equation  (11)  it  is  seen  that  E  is  simply  the  negative  of  the 

n<s  equivalent  vector,  p*cos(T)  -  r*sin(T).  Since  this  is  always 

oriented  parallel  to  the  elevation  drive  axis  there  is  no  need  to 

* 

accelerate  J£  inertially,  but  the  only  requirement  Is  that  E  compen¬ 
sate  for  the  hull  motion.  Torque  is  required  only  to  accelerate 
motor  and  gear  inertia  and  overcome  friction.  The  traverse  axis 
can  be  analyzed  in  a  similar  manner.  The  first  term  in  equation 
(12)  is  q.  This  is  always  oriented  parallel  to  the  traverse  drive 
axis  and  thus  requires  no  torque  to  accelerate  inertially.  But 
the  second  term  in  equation  (12),  tan(E){r*cos(T)  -  p»sin(T)},  is 
oriented  at  a  right  angle  to  the  traverse  drive  axis.  This  is  the 
component  of  hull  motion  which  requires  torque  to  accelerate 
inertially.  It  is  extremely  important  that  these  torque  requirements 
in  terms  of  hull  motion  be  clearly  understood.  For  if  they  are  not,  - 
one  could  end  up  designing  a  control  system  that  would  have  the 


Figure  4 

Equivalent  Base  Motion  Vectors 


function  of  moving  the  hull  with  respect  to  the  sight-line  vector. 

L>t  us  now  proceed  to  the  considerations  of  the  actual  control  system 


CHAPTER  3 


SYSTEM  CONCEPTS 

The  traverse  axis  only  will  be  considered  from  here  on.  This 
can  be  done  because  the  two  axes  are  operated  independently.  All 
the  assumptions  and  conclusions  that  follow  can  be  easily  extended 
to  the  elevation  axis. 

I>causa  of  space  limitation  and  the  extremely  large  traverse 
inertia,  the  size  of  the  notor  that  can  be  used  to  drive  the  traverse 
inert i a  is  very  small,  in  a  maximum  torque  output  sense,  in  comparison 

to  the  size  of  the  traverse  inertia.  This  causes  the  system  to  be 

sluggish  when  compared  to  the  high  speed  servo  mechanisms  that  are 
familiar  to  control  engineers.  This  inherent  sluggishness  is  the 
biggest  obstacle  that  must  be  overcome. 

There  are  two  basic  approaches  to  sight-line  stabilization. 

One  method  is  the  familiar  matching  of  an  output  to  a  command  input. 
The  other  method  matches  the  output  of  an  auxilliary  high  speed 
sight  servo  to  the  command  input,  and  then  uses  the  output  of  the 

sight  servo  as  the  input  to  the  sight-line  stabilization  system. 

The  reason  for  this  will  be  explained  later. 

The  familiar  matching  of  input  to  output  is  known  as  the  dis¬ 
turbed  method.  In  this  method  the  traverse  component  of  the  sight¬ 
line  is  fixed  to  the  traverse  inertia,  (see  fig.  5). 
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Figure  5 

Disturbed  System 
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The  main  control  system  consists  of  the  traverse  load  inertia,  the 
traverse  drive  motor,  the  other  related  components  and  the  necessary 
dynamic  relationships.  Because  the  main  control  system  is  inherently 
sluggish  there  is  a  tendency  for  the  operator  to  apply  manual  correc¬ 
tion  to  help  speed  up  the  system.  This  can  have  a  very  negative 
effect  on  the  overall  performance  of  the  system. 

The  second  method  of  sight-line  stabilization  is  the  directed- 
line  method.  In  this  method  the  sight-line  is  fixed  to  an  auxiliary 
sight  inertia  which  is  very  small  when  compared  to  JT«  Because  the 
sight  inertia  is  small  it  can  be  driven  at  relatively  high  speeds. 

The  command  input  is  applied  to  this  high  speed  sight  servo  and 
matched  to  its  output.  The  output  of  the  sl«r,:it  servo  is  then  used 
as  the  input  to  the  main  control  system,  and  the  output  of  the  main 


control  system  is  matched  to  the  output  of  the  sight  servo  (see  fi 

6). 


figure  6 
Directed  System 


This  has  the  effect  of  matching  the  command  input  to  the  main  control 
system  mtput,  but  it  does  so  in  two  stages.  Because  the  command 
input  is  matched  to  the  output  of  a  relatively  high  speed  system, 
the  tendency  to  apply  manual  correction  is  eliminated. 

There  is  no  conclusive  evidence  as  to  which  of  the  two  pre¬ 
viously  discussed  methods  is  more  effective,  although  the  bulk  of 
the  work,  that  has  been  done  is  on  the  directed  method.  The  attention 


here  will  be  focused  on  the  disturbed  method  because  it  is  the 
simpler  of  the  two  methods,  and  yet  it  still  includes  all  the  basic- 
considerations  of  the  dynamics  necessary  for  a  sight-line 


r 
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t 

stabilization  system.  The  simplicity  also  cuts  down  on  the  compu¬ 
tation  time  that  is  required  tor  the  numerical  optimization  which 
will  be  used  later  on. 


t 
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CHAPTER  4 
DESIGN 

For  a  system  that  is  to  be  optimized,  the  system  need  not  be 
designed  to  meet  certain  performance  specifications  (i.e.  rise 
time,  overshoot,  settling  time,  bandwidth)  as  is  the  usual  case. 

This  is  possible  because  when  the  system  is  optimized  all  the  para¬ 
meters  will  be  adjusted  to  satisfy  the  given  performance  index. 

This  performance  index  may  Include  the  classic  time  response  char¬ 
acteristics  or  nay  be  conq>letely  void  of  them..  Whatever  the  case, 
the  time  response  will  he  shaped  during  the  optimization  along 
with  the  frequency  response.  The  objective  will  therefore  be  to 
design  a  system  that  is  a  reasonable  starting  point  for  the  optimi¬ 
zation  routine.  A  reasonable  system  for  a  starting  point  can  i 
defined  as  a  system  that  does  not  have  anything  critically  wrong 
with  it,  such  as  being  unstable.  The  root  locus  method  of  synthesis 
lends  itself  very  well  to  the  design  of  a  system  with  such  loose 
specifications,  and  is  therefore  what  will  be  used.  It  causes  very 
little  difficulty  to  keep  an  approximate  range  for  a  bandwidth  in 
mind  whwi  using  the  root  loc.us  technique.  The  objective  will  be 
to  design  the  -jystem  for  the  optimization  starting  point  to  have  a 
bandwidth  of  approximately  3C0  radians  par  second.  This  value  was 
obtained  as  a  typical  one  from  the  literature,  however  it  is  by  no 
saar-s  a  value  that  should  be  used  for  all  stabilization  systems. 


A  hydraulic  motor  is  used  to  drive  the  traverse  inertia  becaust- 
rore  torque  can  be  obtained  from  a  hydraulic  motor  than  can  be 
obtained  from  an  electric  motor  of  the  same  physical  si?e  A  partial 
block  diagram  of  a  typical  constant  displacement  piston  t_,pe  hydrau¬ 
lic  motor  and  its  control  valve  is  shown  in  fig  7.  A  list  of  the 
variables  and  a  3et  of  numerical  values  (where  applicable)  is  given 
in  table  1. 

Figure  7 

Fartial  Hlock  Diagram  of  Hydraulic  Motor 


18 


Table  1 

List  of  Motor  Variables 

KP  -  servo  valve  gain  15,000  psi/ma 

Kg  -  pressure  feedback  gain  .002  ma/psi 

GV  -  servo  valve  transfer  function  1/{s^/12002  +  1. 4*3/1200  ->*  1 } 

G,  -  pressure  feedback  compensation  s/{(s/20  +  l)»(s/.5  +1)} 
o 

DM  -  motor  displacement  0.016  in-fyrad 

CP  -  valve  orfice  damping  coefficient  6,200  lb/sec/in^ 

8  -  bulk  modulus  of  fluid  300,000  psi 

-  effective  volume  of  fluid  under  coiqiression  1  in^ 

6  -  actual  motor  angle  in  radians 

ma 

-  motor  torque  in  in-lbs 

-  pressure  applied  to  motor  in  psi 
P,  -  compressibility  flow  in  in^/sec 

Qj3  -  displacement  flow  through  motor  in  ia^/sec 

Qj,  -  total  load  flow  through  metering  orfice  in  in^/sec 

i  -  input  current  to  3ervo  valve  in  ma 

Vc  -  motor  and  valve  configuration  input  current  in  ma 
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C?  is  actually  a  nonlinear  term  which  depends  on  the  flow  rate, 
Q^.  The  value  .given  for  CP  is  for  the  small  signal  region.  The 
valve  input  current  will  saturate  at  some  value,  typically  10  ma. 

The  maximum  supply  pressure  also  hr.s  a  limiting  value  with  3,000 
psi  being  representative  for  the  type  of  system  being  considered. 

The  pressure  saturatic-n  will  be  the  limiting  factor  in  the  perfor¬ 
mance  of  the  system.  The  maximum  pressure  c-f  3,000  psi  and  the 
motor  displacement  of  .016  in  /rad  combine  to  provide  a  maximum 
torque  or  48  in-lb  (*3,000  psi  x  .016  in^)  which  will  be  very  small 
for  the  size  of  load  inertia  being  dealt  with.  A  detailed  develop¬ 
ment  of  a  block  diagram  similar  to  fig.  7  can  be  found  iu  reference 
{8>. 

To  complete  the  block  diagram  of  fig.  7  it  must  be  determined 
ir.  what  manner  the  torque,  *■  ,  drives  the  load.  Consider  the  expanded 
partial  view  of  fig.  3b  in  fig.  3.  The  definitions  of  the  variables 
shown  in  fig.  t  are  given  in  table  2. 
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Table  2 

Definitions  of  Drive  Variables 

q  -  previously  defined  base  yaw  angle 

R  -  radius  of  JT  ring  gear 

r*  -  radius  of  Hj,  pinion  gear 

<p  -  inertial  sight-line  traverse  angle 

6  *  -  motor  angle  needed  to  compensate  for  q 
o 

9  -  component  of  motor  angle  that  is  required  to  inertially 

drive  sight-line  through  angle  ($-q) 

0  -  original  pinion  center 

O'  -  pinion  center  position  due  to  q  input 

O’ *  -  final  inertial  pinion  center 

When  the  base  undergoes  a  yaw  angle  of  q  the  motor  must  acce¬ 
lerate  the  motor  and  gear  inertias  through  an  angle  of  0  ’ ,  and 

m 

overcome  friction  to  keep  eg  oriented  along  the  x  axis  as  is 
required  for  sight-line  stabilization.  Geometrically  this  means 
that  points  A  and  A*  must  remain  coincident.  The  two  variables, 
q  and  6  ’ ,  are  related  by  equation  (13). 

-R.q  =  t*.0  *  (13) 

m 


When  it  is  required  to  move  the  sight-line  from  its  initial 


orientation,  e  ,  to  another  one,  e 
’s’  7  s 


t 


the  motor  must  accelerate 
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Jj,  as  veil  as  the  motor  and  gear  inertias.  must  traverse  an  angle 
<£  while  the  pinion  turns  through  an  angle  6m>  $  and  9^  ars  related 
by  equation  (14).  For  this  case  points  A  and  B  coincide. 

(14) 

When  both  q  and  $  occur  simultane-nisly  need  only  be  driven 
inertiallv  through  an  angle  that  is  the  difference  between  $  and  q. 
Roll  without  clip  conditions  require  that  the  arc  length  from  point 
C  to  point  D  be  the  same  as  the  arc  length  fror>  point  E  to  point 
B',  where  B  and  B*  are  coincident  initially.  This  relationship 
can  be  expressed  by  equation  (15). 

R* (r~q>  “  r- ($+9)  (?5) 

m 

Equation  (15)  can  be  rewritten  in  the  form  of  equation  (16). 

rv6Q  -  (R-iO*$  -  R*q  (16) 

It  is  worth  notine  that  equation  (16)  becomes  equation  (14)  when 
there  Is  no  q  input.  The  gear  ratio,  N,  is  defined  by  equation  (17). 

K  -  R/t*  (17) 

Using  this  relationship  equation  (16)  can  he  rewritten  as: 


r4*  6  =  (R-r)  *  $ 

m 


(18) 


9  =  (N-l)-v  -  N-q  =  N- {  <N-1) -<}>/N  -  q} 

m 

Associated  with  the  motor  shaft  is  a  compliance,  XL.  If  friction 
is  neglected,  the  motor  torque  serves  to  accelerate  the  motor  and 
gear  ineriia,  J  t  and  twist  the  shaft.  The  angle  of  twist  is  the 
difference  between  the  actual  motor  shaft  angle,  6  ,  and  the  angle 

IQd 

that  is  predicted  by  equation  (18) .  The  motor  torque  is  given  by  the 
following  expression. 

t  *  J  -s2-e  t  KL>(9  -  N*({N-1}*$/N  -  q)}  (19) 

c>  m  ma  me 

By  using  the  relationship  in  equation  (19)  and  the  block  diagram  in 
fig.  7,  the  block  diagram  of  fig.  9  can  be  obtained.  Typical  values 
of  the  quantities  introduced  in  fig.  9  are  given  in  table  3  below. 

Table  3 
Load  Variables 

XL  -  motor  shaft  compliance  600  in-lb/rad 

2 

J  -  motor  and  gear  inertia  0.00158  in-lb-sec 

m 

-  traverse  load  inertia  1.44(10*)  in- lb-sec2 


H  -  gear  ratio 

A8  -  motor  shaft  twist  angle 


1000 


The  compressibility  flow,  Q  ,  is  very  small  in  comparison  to  the 
displacement  flow  through  the  motor,  and  can  for  all  practical 
purposes  be  neglected. 

Assuming  q  equal  to  zero,  the  transfer  function  from  to  $ 
can  be  readily  found  following  a  series  of  block  diagram  manipula 
tions  (see  fig.  10),  and  is  given  by  equation  (20). 


Figure  9 


Motor  and  Load 


Figure  10 

Block  Diagram  Manipulation  to  find  <J>/V 
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KP • GV • DM • KI • N • KL 


f ( 1+K? • GV • K. • G, ) • (J  .s2+N.{N-l}-KL-KL2.N-{N-l}*KI)+  1 
*  DO  T 


(20) 


DM2-KI-CP-S.  (^-s^-CN-D-KL) 


I 


where: 


KI 


1/KL _ 

J  -s^/KL  +  1  . 
m 


(21) 


By  substituting  the  appropriate  numerical  values  which  were  given 


previously,  the  transfer  function  of  equation  (22)  is  obtained. 


C 


£ 

V 

c 


_ 9.1?.-(10*2)-(s+20)-(8+.5) _ 

s  •  (s+.  54)  •  (s+4 . 37  )  •  (s+1284 .  58K82+2{ .  21  }370s+3702) . 
(s2+2 • { . 52 } • 1 2C0 • s+12002 ) 


The  last  tarn  in  the  denominator  of  equation  (22)  will  have  a 
flat  frequency  response  out  to  approximately  1200  rad/sec  ,  and  its 
time  response  is  very  fast  when  compared  to  the  other  terms.  For 
all  pratical  purposes  it  can  be  normalized  and  set  equal  to  unity 
since  it  will  have  little  effect  on  either  the  frequency  response  or 
the  time  response.  Dividing  the  numerator  and  denominator  of  equation 
(22)  by  12002  and  assuming  that  the  dyna*iic  effects  of  (s2/12002  + 

2* { .52}*s/1200  +  1)  are  negligible,  equation  (22)  can  be  simplified  to 
equation  (23). 


€ 


A 


* 


f  i;  'I 1  w  ■'  ■  **' ^ 
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±  _ G.34.(106Ms+20).(s+.5) _ 

V  *  s •  (s+. 54).(s+4 . 37Ms+1284.58Msz+2{  . 21  }370s+37(K) 
c 


Consider  the  transfer  function  of  the  servo  valve. 


n\ 1  3  -  i,  ■ 

8^/1200^  +  2. (.7). s/1200  +  1 


(24) 


This  resembles  very  closely  the  term  that  has  just  been  neglected. 

If  GV  is  set  equal  to  unity  when  substituting  the  numerical  values 

into  equation  (20),  the  following  revised  transfer  function  for 

$/V  is  obtained . 
c 

£  _ 6.34«(106)-(8»2nWs+.5) _ 

V  *  s-  (s+.55) *  (s+4.37)*  ( s+1060. .31}406s+4062)  U 
c 


This  ie  in  very  good  agreement  with  equation  (23),  and  therefore 
GV  will  be  assumed  to  be  unity.  There  is  no  way  to  readily 
eliminate  the  term  l/(s  +  1060.5)  in  equation  (25)  by  a  block 
diagram  alteration.  It  must  therefore  be  tolerated. 

Since  $  is  defined  to  be  the  inertial  component  of  traverse 
motion  it  should  be  zero  for  a  yaw,  q,  input  which  is  the  non- 
inertial  component  of  traverse  motion.  This  is  the  same  as  requiring 
that  40  be  zero  for  a  q  input  (see  fig.  9).  Assuming  that  $  does 
equal  zero,  and  setting  equal  to  zero  yields  the  simplified 
block  diagram  in  fig.  11  for  q  inputs  only.  The  validity  of 
assuming  that  <f>  and  equal  zero  will  be  seen  later  on.  G^  and 


Substituting  the  appropriate  numerical  values  into  (31)  yields 
equation  (32). 


G 

<1 


(32) 


liius  if  an  input,  -  N*q,  is  applied  through  a  transfer  function, 

G^,  and  the.  output  is  applied  to  the  stabilization  system  as  shown 
in  fig-  11:  there  will  be  no  motion  of  (•  for  a  yaw  input. 

The  design  of  the  hardware  to  provide  a  transfer  function  for 
0  is  a  difficult  if  not  impossible  task.  The  best  that  can  be 
hoped  for  is  a  good  approximation,  but  it  is  not  the  purpose  of 
this  writing  to  design  hardware.  The  purpose  is  to  show  the  necessary 
mathematical  relationships  for  sight-line  stabilization. 

Yaw  Inputs  have  a  negligible  effect  on  the  performance  of  a 
sight-line  stabilization  system  when  compared  to  the  inertia?,  inputs. 
Therefore  further  efforts  will  be  concentrated  on  the  inertial  inputs 
only. 

Consider  again  fig.  9  with  the  two  q  inputs  and  the  compress¬ 
ibility  term  removed,  and  GV  set  equal  to  unity.  It  is  repeated 
below  as  fig.  12  for  convenience.  The  transfer  function  for  fig.  12 
was  given  by  equation  (25)  which  is  also  repeated  below  as  equation 
(33). 


There  are  a  number  of  rate  gyros  available  that  have  a  natural 
frequency  much  larger  than  300  rad/sec.  For  a  rate  gyre  of  this  type, 
the  output  will  ca  approximately  proportional  to  the  rate  only. 


X’S 


ma-sec/raa 


(35) 


Similarly  for  a  posit?  >n  gyro  with  the  same  restrictions  on  natural 
frequency,  its  transfer  function  can  be  approximated  by  a  simple 
gain. 


£ 

»  — tttt  n~*7T  y  -  K  me/rad  (36) 

Pv.  f  2?sAc  +  I 

It  must  cow  be  determined  whar  KV  and  should  be.  A  pole  zero 
plot  for  equation  (33)  is  given  in  fig.  14.  The  scale  has  been  ex¬ 
panded  in  the  region  ;;ear  the  origin  in  order  that  the  smaller  terms 


can  be  shown. 


Figure  14 

Pole  ?.ero  Plot  of  Equation  (33) 


domewhat  arbitrarily  it  wi.il  be  assumed  that  KV  is  unity.  It  should 
be  remembered  that  this  is  only  a  starting  value  for  optimization. 

It  is  easy  to  see  fro®  fig.  1*.  that  the  five  poles  and  zeroes  near 
the  origin  will  dominate  the  time  response.  When  plotting  the 
root  locus  for  fig.  13  the  velocity  feedback  will  cancel  out  the 
pole  at  the  origin,  and  the  pole  at  0.55  can  be  assruoed  to  cancel 
with  the  zero  at  0.5.  It  is  reasonable  to  assume  that  G^c  contain 
a  term  (s  +  4.37)/{s  +  20)  for  cancellation.  If  this  and  a  gain 
were  all  that  were  included  in  G„_,  then  the  root  locus  would  tend 


co  the  right  half  plane,  and  there  would  be  little  room  for  the 
bandwidth  of  the  system  to  be  increased  without  lowering  the  damping 
ratio  of  the  complex  pole  pair.  3y  introducing  two  lead  networks 
of  the  form  (s  +  150) /(s  +  2000)  the  root  locus  will  take  the  form 
of  fig.  15. 


Figure  15 


Root  Locus  of  Velocity  Loop 


Tiie  first  thing  to  note  is  that  the  two  zeros  at  150  are  too  small 
if  the  frequency  response  is  to  be  flat  out  to  approximately  300 
rad/sec.  Also  a  closed  loop  pole  will  fall  someplace  between  0  and 
150.  Therefore  the  two  zeroes  must  be  increased  in  -magnitude  to 
help  draw  the  root  locus  to  the  left.  After  some  trial  and  error, 

700  was  chosen  to  be  n  suitable  value  for  the  zeroes.  Thus  G  will 

iv 

contain  two  lag  networks  of  the  form  (s  -i  70C)/(s  +  ISO).  If  these 
two  lag  networks  and  a  gain,  PK,  were  all  that  were  contained  in  Gpc> 
the  root  locus  of  fig.  16  would  tend  to  the  right  half  plane  from 
the  smaller  complex  pole  pair,  and  there  would  be  little  room  fcr 
the  damping  ratio  to  decrease  during  optimization,  if  desired.  This 
can  be  averted  by  introducing  a  lead  network  such  as  (e  +  700)/ 

(s  +  3J00)  into  Gpc>  A  root  locus  for  the  position  loop  with  the 
previously  discussed  compensation  is  given  in  fig.  18.  For  a  position 
compensator  gain,  PK,  of  1200,  $/R  will  be  given  by  equation  (39). 


± 

R 


_ 6.23-(1012)-(s  +  7QO;3 _ 

(s  +  297) -(s  +  145  +  j  *540) • (s  +  145  -  j-540)-' 
(a  +  3303  +  j *584)* (s  +  3303  -  j-584)- 
(s  +  559  +  j’1314)*(s  +  559  -  j-1314) 


(39) 


Therefore  it  can  be  assumed  that  a  reasonable  starting  value  of  G 

j  w 

is  given  by  equation  (40) . 


PC 


120C- 


(s  *  7QQ) J 


<3  +  150)2, (s  +  3000) 


(4G) 
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The  transfer  function  $/R  satisfies  the  300  rad/sec  bandwidth 
requirement ,  and  the  smallest  damping  ratio  present  is  the  0.26 
associated  with  the  complex  pair  (s  +  145  +  j*540)*(s  +  145  -  j • 

540).  Although  a  damping  ratio  of  0.26  gives  approximately  45" 
overshoot  the  system  will  have  nowhere  near  this  amount  of  overshoot 
due  to  the  pressure  saturation  previously  discussed,  nor  will  the 
system  be  anywhere  near  as  fast  as  the  denominator  of  equation  (39) 
might  indicate.  The  whole  sight-line  stabilization  system  is  shown 
in  fig.  19.  It  includes  inertial  inputs,  R,  as  well  as  the  non- 
inertxal  input  q.  This  is  not  the  final  system,  but  only  a  starting 
point  for  the  optimization. 

It  can  be  seen  that  for  a  q  input  only,  R  will  be  equal  to  zero, 
if  depends  only  on  what  is  between  It  and  the  input  to  the  system, 
in  this  case  the  only  input  is  q.  Since  was  designed  to  cause 
if  to  be  zero  for  q  inputs,  there  will  be  no  position  feedback  from 
if.  Therefore  will  equal  zero,  and  it  is  seen  that  the  assump¬ 
tions  made  earlier  in  regard  to  and  q,  when  designing  G^,  were 
valid. 

Although  a  strict  linear  analysis  was  used  throughout  this 
chapter,  and  there  was  no  non-linear  check  on  the  stability  of  the 
system  due  to  the  pressure  saturation;  it  is  belived  that  there 
will  be  no  problem  with  stability.  This  is  because  of  the  limiting 
value  on  the  torque  of  the  system. 


W . trtjfcMjkl, 


Final  System 
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CHAPTER  5 

SIMPLIFIED  VERSION 

Numerical  computation  time  can  be  a  prohibitive  factor  when 
optimizing  a  system,  as  is  the  case  for  the  system  of  chapter  4.  It 
is  the  objective  of  this  chapter  to  develop  a  simplified  model  that 
will  reduce  the  amount  of  computation  time  required  for  optimization. 

An  examination  of  fig.  19  reveals  an  undamped  second  order  term 
in  the  expression  for  the  motor  load.  This  expression  is  6iven  below. 


1/KL 


J  -s2/KL  +  1 

Cl 


(4n 


Multiplying  the  numerator  and  denominator  of  this  expression  by 
KL/J^,  and  substituting  the  numerical  values  given  yields  (42). 


632.9 


s2  +  616. 22 


(42) 


The  numerical  optimization  routine  that  will  be  used  later 
requires  that  the  differential  equations  of  the  system  be  integrated 
repeatedly.  When  numerically  integrating  an  expression  with  a 
natural  frequency  as  high  as  the  one  in  (42)  a  large  amount  of 
compute  time  is  required  because  the  integration  interval  has  to  be 
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s 
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A 
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made  very  small  in  order  to  achieve  reasonable  accuracy.  For  all 

2 

practical  purposes  the  term  J  *s  /KL  can  be  assumed  to  be  negligible, 
and  expression  (41)  can  be  approximated  by  1/KL.  This  alleviates 
the  need  to  integrate  expression  (41),  however  there  are  other 
expressions  that  also  must  be  dealt  with  in  order  to  reduce  compu¬ 
tation  time  for  the  whole  system. 

Again  referring  to  fig.  19  if  the  loop  containing  the  terms 
2 

N,  N-l  and  N-KL/ (J^* s  )  is  collapsed,  the  following  expression  can 
be  readily  obtained: 

N-KL/J., 

s2  +  h-(n-d*kl/jt  (43) 

By  substituting  the  appropriate  numerical  values  into  (43) ,  expression 
(44)  is  obtained. 


41.66 

2  2 
s  +  209 


(44) 


This  term  also  requires  a  very  small  integration  interval  in  order 
to  achieve  reasonable  accuracy. 

Expression  (43)  arises  out  of  the  dynamical  considerations  for 
motor  shaft  compliance.  If  it  is  assumed  that  the  motor  shaft  is 
rigid,  then  it  is  easy  to  see  that: 


$ 


6 

m 


N-l 


(W>) 
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Neglecting  friction,  the  motor  torque  serves  to  accelerste  the  motor 
and  gear  inertia  through  the  angle  6^,  and  the  traverse  inertia 
(divided  by  the  gear  ratio)  through  ar:  angle  This  relation  is 
expressed  in  equation  (46) . 


m 


T  2  n 
J  *s  *0 
m  m 


+  J  •s*'»4i/N 


(46) 


If  equation  (45)  is  substituted  into  equatxon  (46) ,  and  K-.l  is 
assumed  to  be  approximated  by  *1,  equation  (47)  is  obtained. 
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(47) 


Combining  equation  (47)  with  the  block  diagram  in  f:.g.  7  and  neglect¬ 
ing  compressibility  flew  and  GV  as  was  done  earlier,  the  block  diagram 
in  fig.  20  is  obtained  where: 


J.  -  J_/N  +  J  *N  (48) 

LI  m 


By  ’.sing  the  block  diagram  In  fig.  20  and  substituting  in  the 
proper  numerical  values,  the  expression  for  $/V r  is  obtained^ 


Figure  20 

Simplified  Version  of  Motor  Load  Relation 


Following  the  same  procedure  as  was  used  for  the  development  of 
the  model  in  chapter  4,  the  next  step  is  to  add  velocity  feedback. 
This  is  represented  in  fig,  21. 


Figure  21 
Velocity  Feedback 
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Let  K.  be  unity-  ana  let  G  be  used  for  the  cancellation  of  the 
v  '  v 

seller  terms  in  equation  (49).  is  then  expressed  as  (50). 

VK* {n  ±  .55) -(a  *  4.38) 

Y  *  (s  +  .5)*  (s  +  20) 


VR  is  a  variable  gain  to  be  determined.  'Hie  pole  at  .55  anu  the 
zero  at  .5  were  assuoed  not  to  cancel  because  the  pole  at  .55  was 
fou^d  to  be  somewhat  sensitive  to  changes  in  the  other  parameters 
when  some  preliminary  programs  of  thf.  system  were  run.  By  using 
equation  (50)  the  transfer  function  for  $/V  is  fo>md  to  be; 


i  _  _ VK-L5.02 _ 

V  “  s(s  +  414.90  +  VK*  15.02)“ 
v 


Consider  $/V  with  position  feedback  as  is  rhovm  in  fig.  22. 


Figure  22 


Position  Feedback 


rt 
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Lit  Cy  be  simply  &.  gain,  PK.  Then  the  transfer  function  for  fig. 22 
ia  given  by  equation  (52). 


i. 

R 


PK-VIM5.02 


<5  +  (414.90  +  VK'15.02)  •  «  -*• 


(52) 


The  expression  414.9  +  YX*15.02  will  correspond  to  the  2cm  term  of 
a  quadratic.  Remembering  that  a  bandwidth  of  approximately  300  is 
desired,  u  can  be  chosen  so  as  to  satisfy  this  requirement.  Let 
u  be  315,  also  let  £  be  unity  tor  an  initial  value.  By  substituting 
these  values  in  equation  (52),  VK  is  found  to  be  14.32  and  PK  is 
found  to  be  461.29. 

Fig  23  is  then  the  final  simplified  configuration  for  a  system 
to  be  used  as  a  starting  point  for  optimization.  This  system  is 
for  inertial  inputs  r-i-jy,  i.e.  R  is  the  inertial  component  only.  The 
yaw  inputs,  q,  v^re  ignored  altogether  because  in  this  simplified 
version  they  are  trivial. 
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A  recsonable  question  to  ask  is:  How  good  is  the  approxima¬ 
tion  of  the  detailed  model  that  is  given  by  tha  simplified  model 
in  chapter  5?  The  answer  to  this  question  can  be  found  by  running 
the  model  of  chapter  a  with  a  very  large  number  inserted  for  KL. 

This  approximates  the  case  of  a  rigid  motor  shaft.  However  if 
equation  («3)  is  considered  it  can  be  seen  that  as  KL  is  increased 
so  is  the  natural  frequency  of  the  undamped  second  order  term, 
and  this  causes  the  computation  time  for  the  solution  of  the  system 
equations  to  go  up.  An  attempt  was  made  to  run  the  system  with  a 

g 

very  large  value  for  KL(10  ).  but  the  amount  of  computation  time 
required  for  this  was  prohibitive.  However,  the  simplified  system 
is  thought  to  give  a  good  approximation  to  the  detailed  version 
because  the  transient,  response  of  the  simplified  system  (see  appendix 
A  and  Fig.  26;  is  on  the  order  of  the  transient  responses  er countered 


in  tho  literature 


CHAPTER  6 


OPTIMIZATION 

The  term  optimization  is  an  abstract  one  to  say  the  least. 

There  are  many  different  definitions  of  optimization.  Basically 
to  optimize  a  system  means  to  minimize  a  performance  index  (cost 
index)  while  the  system  goes  from  one  state  to  another.  This  per¬ 
formance  index  is  a  function  of  the  parameters  of  the  system.  The 
system  parameters  include  all  the  physical  variables  of  the  system 
(gains,  compensator  values,  inertias,  etc.);  and,  in  the  strictest 
sense,  the  input.  The  input  will  not  be  assumed  to  be  a  parameter 
here.  This  corresponds  to  an  input  over  which  there  is  no  control, 
a  case  that  is  not  uncommon  in  the  control  system  area. 

There  are  many  forms  for  the  performance  index.  Integral  forms 
are  a  very  popular  type.  Among  the  integral  fcrme  is  the  integral 
squared  error  (ISE).  If  a  system  goes  from  one  state  at  time  t^ 
to  another  state  at  time  then  the  ISE  is  given  by  equation  (53), 
where  R  is  the  Input  to  the  system  and  C  is  the  output  of  the  system. 


ISE 


(R-C)2dt 


(53) 


With  the  use  of  Parceval’s  theorem  {7}  and  tables  {7},  the  ISE  can 
be  computed  with  some  ease  when  *»  0  and  t^  *  ®. 


This  is  not  to 


say  that  a  minimum  ISE  is  easily  obtained  analytically. 

Numerical  optimization  routines  have  gained  in  popularity 
because  of  the  difficulty  in  applying  analytical  techniques  in  opti¬ 
mization.  The  key  to  a  numerical  technique  is  the  method  used  to 
vairy  the  system  parameters  in  the  search  for  the  optimum  solution. 
There  are  two  classes  of  search  techniques:  indirect  search  and 
direct  search.  All  indirect  search  methods  employ  the  use  of  the 
gradient  of  the  performance  surface  in  one  manner  or  another. 

Highly  oscillatory  or  even  discontinous  performance  surfaces  are 
very  common  in  the  control  system  area.  This  highly  oscillatory 
nature  can  cause  a  gradient  to  be  calculated  incorrectly ,  thus 
affecting  the  optimization.  For  this  reason  direct  search  methods 
are  more  suited  to  optimal  control  prubl-ims. 

There  are  many  fine  direct  search  techniques.  Among  them 
are  the  methods  of  Hook  ard  Jeeves  {4}  and  the  method  of  Rosenbrock 
{9}.  Lange-Nielsen  {6}  has  modified  Rosenbrock' s  method  and  used 
it  to  develop  an  optimization  program.  A  detailed  description  of 
Rosenb rock’s  method  and  Lange-Nielsen' s  modifications  can  be  found 
in  reference  [63.  A  step  by  step  implementation  of  Lange-Nielsen ’ s 
program  can  also  be  found  in  reference  [6], 

Lange-Nielsen 's  program  will  be  used  to  optimize  the  sight-line 
stabilization  system  of  chapter  5  using  ISE  as  the  performance  index. 
In  Lange-Nielsen ' s  program  the  system  parameters  are  perturbed  using 
the  modified  version  of  Roscnbrock's  method,  and  after  each  pertuba- 
tion  the  system  differential  equations  are  solved  using  CSMP  {5}. 
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The  performance  Index  is  calculated,  and  the  procedure  Is  repeated 
until  convergence  is  cssured.  The  reader  is  referred  to  reference 
CiOj  for  a  verification  of  the  accuracy  of  Lange-Nielsen *s  program. 

Before  the  actual  optimization  can  take  place,  the  parameters 
that  are  to  be  varied  must  be  decided  upon.  The  procedure  that  will 
be  used  here  will  be  one  of  elimination  of  inappropriate  parameters. 

Consider  figure  23.  Let  the  zeroes  of  the  compensator  be  known 
as  ZL1  and  ZL2,  and  let  the  poles  be  known  as  PL1  and  PL2.  Also  let 
the  poles  of  G -  be  called  PI  and  P2,  All  the  system  oarameters  can 

U 

now  be  listed  in  table  4. 


Table  4 


List  of  System  Parameters 


PK 

PL2 

DM 

VK 

KP 

CP 

ZL1 

K6 

N 

ZI.2 

PI 

JT 

Pll 

P2 

J 

m 


It  Is  intuitively  obvious  that  Jffi  and  will  tend  to  their 
lower  limits,  and  DM  will  tend  to  its  upper  limit.  If  the  expression 
for  JT  is  considered,  it  can  be  seen  tiiat  N  will  tend  to  a  very 
large  number  in  order  to  minimize  J .  It  can  therefore  be  assumed 

w 

that  the  bshavior  of  these  four  parameters  during  optimization  is 
known,  and  that  the  values  given  have  already  been  extended  to  their 
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limits.  CP  is  a  small  signal  value  for  a  nonlinear  term,  and  there¬ 
fore  to  use  it  as  a  parameter  would  have  little  meaning.  PI  and  P2 
are  part  of  the  pressure  feedback  transfer  function,  G^.  It  is  not 
always  possible  to  order  hardware  that  has  the  exact  physical 
characteristics  that  are  wanted.  This  may  be  the  case  with  PI  and 
P2,  therefore  they  will  also  be  deleted  as  parameters.  PL1  and  PL2 
are  see*  lo  cancel  PI  and  P2  directly,  and  since  PI  and  P2  are  de¬ 
leted  as  parameters,  PL1  and  PL2  will  also  be  deleted.  KP  and 
are  part  of  the  servo  valve  pressure  feedback  configuration.  But 
since  they  are  gains  it  will  be  assumed  that  they  can  be  altered 
without  much  trouble.  Some  preliminary  runs  of  the  program  have 
shown  that  KV  remains  approximately  equal  to  unity.  KV  will  there¬ 
fore  be  left  fixed  at  one.  There  are  no  foreseen  reasons  for  not 
using  the  remaining  variables,  therefore  the  variables  to  be  used 

in  the  optimization  are  KP,  PK,  VK,  K, ,  ZL1  and  ZL2.  In  order  that 

o 

the  compensator  be  put  in  proper  fora  for  CSMP,  VK  must  be  modified 
as  follows: 


VK  *  14.32-PL1-PL2/ZL1/ZL2  (54) 

The  parameters  are  listed  below  in  table  5  along  with  their  start¬ 
ing  numerical  values  for  convenience. 

Everything  is  now  ready  for  the  actual  optimization,  which 
is  done  with  the  use  of  the  optimization  program  developed  by  T. 
Lange-Nielstn  [10]as  previously  mentioned.  The  optimal  parameters 
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are  listed  belov  in  table  6  for  several  magnitudes  of  step  inputs 
along  with  the  percent  overshoot  of  the  system  and  the  optical 
value  of  ISE.  The  step  magnitudes  are  representative  of  those  en¬ 
countered  in  pratice.  The  changes  In  percent  overshoot  are  due  to 
the  non-linearity.  A  sensitivity  eheck.  of  the  system  with  the  0.1 
step  is  given  in  t~bie  7. 


Table  5 


Initial  Values  of  Parameters  Used  in  Optimization 


KP  =  15000  psi/na 


0.002  ma/psi 


PK  *  461.29 


ZL1  «  0.55 


VK  -  3.43 


ZL2  «  4.36 


Table  6 


Optimal  Parameters 


Step  Magnitude  |  0.05  {  0.10  0.15 


29231.1 


29170.8 


.0018154  .0025239  .0006002  i .0033426 


540.687  463.931  435.085  ! 400.212 


3.60303  4.14173  2.93995  j  2.9 


.548911  .512743  .497890  1.547466 


9.05973  6.06801  5.74567  14.39803 


28.2  15.87 


2. 4763E-4  1. 3941E-3  3. 8416E-3  I 8.0115E-3 


+10 


1  3945E-3 


. 

* 

KP 

l 

-10 

1.3941E-3 

5.59 

+20 

1.3945E-3 

7.48 

O 

<N 

1 

1.3945E-3 

4.22 

K6  ' 

+10 

1.3941F.-3 

5.59 

-10 

1.3945E-3 

7.48 

+20 

1.3946E-3 

4.05 

-20 

1 

1.3945E-3 

7.48 

• 

PK 

+10 

1.3945E-3 

7.48 

-10 

1.3948E-3 

3. 

73 

+20 

1.3945L-3 

7.48 

-20 

1.3948E-3 

3.73 

+10 

1, 3945E-3 

7.48 

wy 

-10 

1.3945E-3 

4.14 

YJc 

+20  . 

1.3.'45E-3 

7.48 

-20  i 

1.3948E-3 

3.; 

73 

HO 

1.4043E-3 

13.: 

74 

zu 

-10 

1.4041E-3 

+20 

1.4285E-3 

19.15 

-20 

1.4314E-3 

— 

Parameter  I 


Table  7 
(Continued) 


Z  Perturbation 

ISE 

+10 

I.3950E-3 

-10 

1.3951E-3 

+20 

1.3957E-3 

-20 

1.3959E-3 

The  frequency  response  of  the  system  that  had  the  0.1  step  as  the 
system  input  is  found  in  figures  24  and  25 ,  and  the  transient  response 
of  the  same  system  Is  found  In  figure  26  and  appendix  A .  A  discussion 
of  all  the  results  obtained  in  this  chapter  and  the  conclusions  that 
can  be  dram  from  these  results  is  found  in  chapter  7. 


CHAPTER  7 


DISCUSSION  AND  CONCLUSIONS 


The  performance  of  the  system  fellows  an  expected  trend  when 
compared  to  the  trend  of  the  inputs.  It  can  be  seen  from  table  6 
that  as  the  magnitude  of  the  step  input  decreases,  the  ey3tem  over¬ 
shoot  Increases.  Ac  the  magnitude  of  the  step  decreases  there  is 
less  torque  requited  to  drive  the  aye ten;  hence  there  is  less  satur¬ 
ation  of  the  load  pressure,  and  the  system  behaves  more  like  a  lin¬ 
ear  system. 

The  important  question  is:  Which  set  of  parameter  values  should 
be  used  for  the  final  system?  Since  the  physical  inputs  co  the  system 
will  be  of  a  random  nature  (pitch,  roll  and  yaw),  there  will  be  a 
wide  range  of  inputs  to  the  system.  The  system  will  obviously  not  be 
able  to  perform  in  an  optimal  meaner  for  all.  the  inputs.  One  approach 
to  determine  the  beet  eet  of  parameters  might  be  to  use  each  step  as 
the  input  to  each  optimal  system  and  compare  the  results. 

A  consideration  of  the  sensitivity  check  in  table  7,  and  an 
overall  view  of  the  parameters  in  table  6  may  serve  just  as  well. 

If  the  change  in  the  percent  overshoot  from  the  optimum  value 
(.-.ptimum  value  i*  5.87  X)  is  considered,  then  the  system  is  some¬ 
what  insensitive  to  all  the  parameters  except  for  ZL1.  An  examination 
of  the  values  for  ZLl  given  in  table  6  shows  that  the  smallest  value 
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for  ZLl  is  approximately  91  %  of  the  largest  value.  If  a  value  some¬ 
where  in  between  is  chosen,  the  performance  of  the  system  for  any 
of  the  inputs  will  not  be  affected  greatly. 

EF,  K,.  and  ZL2  all  have  wide  variations  in  their  parameter  values 
o 

as  can  be  seen  from  table  6,  however  an  examination  of  table  7  shows 
thee  these  three  parameters  are  all  very  insensitive  with  respect  to 
the  change  in  percent  overshoot.  PK  and  VK  do  not  vary  a  great  deal, 
nor  are  they  very  sensitive. 

Considering  the  data  in  chapter  5,  the  final  decision  on  the 
parameter  values  must  be  made  only  after  considering  all  the  available 
facts.  Facts  such  as  linearity  of  elements  (amplifier  gains,  compen¬ 
sator  poles,  etc.),  and  the  actual  shape  of  the  input  distribution. 
Ideally,  if  the  actual  system  input  were  available,  the  optimal  para¬ 
meters  could  be  found.  However  the  actual  system  input  la  not.  always 
available  as  is  the  case  here.  Therefore  the  optimal  parameters  of 
chapter  6  can  only  be  used  as  one  factor  in  the  final  choice  for  tne 
final  system. 

The  frequency  response  plots  of  figures  24  and  25  reveal  that 
;he  30C  rad/sec  bandwidth  requirement  has  been  satisfied  for  the 
optimal  system  with  the  0.1  step  input.  They  also  reveal  that  there 
is  no  problem  with  stability  as  is  sometimes  encountered  in  systems 
with  non-linear  elements.  The  transient  response  in  appendix  A  S  Fig.  26 
proves  the  dominance  of  the  pressure  saturation  over  the  system  per¬ 
formance  as  was  expected.  One  final  conclusion  can  be  drawn  by  noting 
the  si^rJ.1  amount  of  overshoot  for  the  systems  with  the  three  l_rgest 
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inputs,  ana  concluding  that  for  systems  with  strict  lltiitstrcns  on 
available  drive  torque  the  overshoot  tends  to  be  small. 


Append!:?  A 


System  Output 

A  list  of  the  print-plotted  variables  and  their  definitions 
3re  given  in  table  8. 

Table  8 

Output  Variables 

PHI  -  system  output 
X9Q0  -  torque 

ISE  -  integral  squared  error 
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